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Abstract 

Photon energy moments in B — > X s + 7 and the impact of experimental cuts are analyzed, 
including the biases exponential in the effective hardness missed in the conventional OPE. 
We incorporate the perturbative corrections fully implementing the Wilsonian momen- 
tum separation ab initio. This renders perturbative effects numerically suppressed while 
leaving heavy quark parameters and the corresponding light-cone distribution function 
well defined and preserving their physical properties. The moments of the distribution 
function are given by the heavy quark expectation values of which many have been ex- 
tracted from the B-^X c iv decays. The quantitative estimates for the biases in the heavy 
quark parameters determined from the photon moments show they cannot be neglected 
for E cut > 1.85 GeV, and grow out of theory control for E cut above 2.1 GeV. Implications 
for the moments in the B~^X c lv decays at high cuts are briefly addressed. 
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1 Introduction and Basics 



The decays of beauty hadrons - recognized for a long time as an area allowing high 
sensitivity probes of fundamental dynamics - have entered now also the realm of precision 
physics on both the experimental and theoretical side. Inclusive semileptonic decays of B 
mesons allow the accurate extraction of the CKM parameters and V u b through their 
integrated rates and of other basic heavy quark parameters through the moments of their 
spectra. These heavy quark parameters include the heavy quark masses m& and m c which 
are 'external' to QCD, and B expectation values of local heavy quark operators like the 
kinetic /4 or chromomagnetic /j,q. The latter are controlled by the strong forces and thus 
can in principle be calculated within QCD once theoretical control has been established 
over nonperturbative dynamics. The motivation for knowing the heavy quark parameters 
accurately goes beyond wanting to extract V c b and V u b- their values provide challenging 
tests for the available treatment of strong dynamics. There has been significant recent 
progress in this direction both on the theory and experiment [1, 2, 3] side, with BABAR 
establishing new benchmarks by combining the high statistics and quality of their data 
with a robust theoretical treatment. 

Detailed studies of the shape of the photon spectrum in B^>X S +^ transitions produce 
complementary information: (i) They allow a systematically independent determination 
of fi^ and maybe even higher order terms, (ii) They lead to a more direct and thus po- 
tentially more accurate extraction of m& than in semileptonic distributions which depend 
on a combination of m& and m c . (iii) They provide immediate access to the distribution 
function of the b quark inside the B meson, the knowledge of which reduces the model 
dependence in extracting V uh from inclusive B — > X u tv transitions. 

Claims have recently appeared in the literature challenging the classical OPE results 
relating the moments of the inclusive distributions to the universal local heavy quark 
expectation values in the decaying hadrons. They assert that the expectation values 
extracted from the b-^civ decays cannot be related to the similar moments of the light- 
cone distribution function controlling the heavy-to-light decays like b — > s + 7, due to 
intervention of the perturbative effects. We find no justification for such revisions, and 
contend that the perturbative corrections to the moments in both types of the decays 
can be treated in the standard OPE approach, see Ref. [4] for a recent analysis. The 
heavy quark parameters extracted from B^>X c tv definitely can be used to constrain the 
light-cone distribution function. 

In practice, when measuring the photon energy spectrum in B^>X S + ^ experimenters 
have to impose lower cuts on the photon energy to suppress backgrounds. Those cuts are 
presently in the range 1.8 to 2.1 GeV with little realistic hope for lowering them in the near 
future. Yet, as emphasized recently [5], such cuts degrade the treatment of inclusive B 
decays in the short-distance expansion. High cuts give rise to a new mass scale parameter 
Q in the expansion, which we refer to as 'hardness', along with the original parameter 
m b . Say, for the inclusive B^>X S +^ decays we have 



Q^M B - 2E t 



'cut 
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Some nonperturbative contributions are then described in inverse powers of Q rather than 
m&. For high E cnt Q becomes too small for a 1/Q expansion to be reliable. 

The usually employed practical version of the operator product expansion (OPE) 1 
often manifests this deterioration in the power corrections explicitly - in particular for 
semileptonic decays - through the Wilson coefficients of the higher-dimensional operators: 
they are functions of E cut /rrib, which effectively transform l/m\ into l/Q k terms. A similar 
reduction of the effective hard scale is present also in the purely perturbative contributions 
and can be traced through a decrease in the effective normalization scale for a s , once the 
BLM corrections are included. 

The OPE-based theory usually leaves out, however, exponential terms like e~ M ^ had 
with /^had denoting the typical momentum scale of nonperturbative QCD, and M stands 
for the generic high mass scale used in the expansion. Those are insignificant as long 
as M is driven by the b quark mass scale and exceeds a few GeV. Yet, as the presence 
of high cuts causes Q to decrease towards 1 GeV, the terms which are exponential in Q 
rather than in m& could quickly become significant. 

The simplicity of two-body kinematics for s + 7 prevents the cuts from affecting 
the Wilson coefficients of the higher-dimensional operators. They do not directly reveal, 
therefore deterioration of the expansion parameter with high cuts. The cut dependence 
will surface and manifest this once the perturbative corrections to these power terms 
are incorporated, yet this still is seen as a future level for theory. Let us note that 
these exponential effects are not connected with the potential growth of the perturbative 
corrections to the Wilson coefficients when the cuts are raised. We have found indications 
that in the actual B decays such perturbative corrections are insignificant, once a physical 
renormalization scheme is used. 

The pilot study of Ref. [5] has shown that using the conventional OPE expressions 
which lack such exponential terms introduces substantial biases into the values of m b and 

extracted from truncated photon energy moments - i.e. moments evaluated with a 
cut on the photon energy - for E cnt ~ 2 GeV. Furthermore, applying the corresponding 
bias corrections, even estimated in a simplified manner, resulted in a surprisingly good 
agreement in the OPE description between 6^5 + 7 and 6— >c + £v, and thus removed an 
otherwise apparent inconsistency. Therefore, in the present paper we want to present a 
more careful theoretical evaluation of the photon moments with cuts, even if it is somewhat 
model-dependent, and to put them into a form suitable for a global analysis of various B 
decay data. 

It had been appreciated quite some time ago that the usual application of the OPE to 
radiative B decays must be incomplete, since it predicts the nonperturbative component of 
the photon energy moments to be independent of E cnt . Ref. [7] tried to assess the potential 
effect of the biases, in our terminology, by applying general inequalities based only on the 
positivity of the distribution function. Being very general, the inequalities however are 
saturated by quite unphysical functions having a support consisting of isolated points. 
Therefore they are too weak to be of interest in practice. Furthermore, these inequalities 

1 A detailed discussion of the principal simplifications employed in what we conventionally refer to as 
the 'practical OPE', can be found, e.g. in the review by M. Shifman [6], and references therein. 
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do not generally hold in the presence of the perturbative corrections. Hence this approach 
has limited utility. 

In Ref. [7] it was also stated that the inequalities are violated once the cut is placed too 
high, and this was interpreted as physical evidence for the OPE breaking down for high 
cuts. We do not think such an interpretation can be correct. The discrepancy noted there 
has a purely algebraic nature. It was actually caused by the Darwin term in Eq. (5.11) of 
Ref. [7] entering with the wrong sign. Correcting it leads to the inequality being trivially 
satisfied at arbitrary E cut and thus eliminates the contradiction, as it has to be for a 
purely arithmetic bound. 2 

Alternative numerical estimates based on the AC 2 M 2 ansatz [8] led to a gross un- 
derestimate: the biases emerged as negligible at cuts below 2 GeV. The perturbative 
corrections were not incorporated at this point, and too low nonperturbative parameters 
were used as a result of following the pole scheme of HQET. 

Similar complications plagued attempts to provide definite numerical estimates in 
Ref. [9]. The analysis there did not use consistently defined Wilsonian heavy quark 
parameters, and operated in terms of somewhat indefinite HQET analogies. The pertur- 
bative corrections in that scheme are not stable and significantly change subtle effects 
like biases in the photon moments. On the other hand, the values of the nonperturba- 
tive parameters determining the assumed distribution function were varied in an overly 
wide range including unphysically small values. Therefore it was difficult to arrive at 
conclusive estimates, and Ref. [9] came up with no certain result: the estimates ranged 
from the biases being totally negligible even at E cut = 2 GeV to very significant down to 
£ cut = 1.8GeV. 

As described in Ref. [5], relying on the Wilsonian implementation of the OPE allows 
to derive more robust statements. The perturbative corrections to physical effects are 
moderate, and the nonperturbative parameters have definite values, often known with 
small numerical uncertainties. For instance, BaBar has recently presented the results of a 
comprehensive fit to the moments they have measured in inclusive semileptonic B decays, 
obtaining [2]: 

m 6 (lGeV) = 4.61 ±0.05|cxp±0.04| H QE±0.02| th GeV 

/4(lGeV) = 0.447 ±0.035| cxp ±0.038| HQ E±0.010| th GeV 2 , (2) 

even without incorporating independent theoretical constraints on the heavy quark pa- 
rameters (see, e.g., [10]). The results are in good agreement with earlier preliminary 
DELPHI values [1, 11]. 

The analysis of B — > X s + j is complicated by the fact that it is induced not only by 
the genuinely local weak operator 

7 = ^m b sa, u F^(l+ l5 )b, (3) 



2 The bound follows from the positivity of the integrals of F(z) (1— z)(l— tz) 2 for arbitrary t, which is a 
quadratic polynomial in t. Its determinant is therefore negative or zero, which yields the inequality Since 
the difference 1 — P n (E cu t), with P n the weight used in Ref. [7] belongs to the above class of functions 
(t=l), Eq. (5.11) is trivially satisfied. 
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but depends also on additional terms in the effective Lagrangian 

AC 8 

^ = ^E^ ( 4 ) 

(we use the standard definition of Oj which is given, e.g., in Ref. [9]). In particular, decays 
mediated by the operator O2 = |s fc 7 M (l — 75)6' c'7 M (l — 7s)c fe (with the superscripts k, I 
denoting color indices to be summed over) contribute through a (real or virtual) cc-pair 
converting into a photon. The effect of operators other than O7, most notably of 2 , 
is essential for a precise evaluation of the decay rate [12]. However, it is expected that 
the normalized photon energy moments are not changed noticeably. The corresponding 
refinements of Refs. [13, 9] modify them by negligible amounts and therefore are irrelevant 
in practice. While illustrating their effect for completeness, we need not dwell on the 
theoretical subtleties associated with their evaluation. 

The effects not captured directly by the standard OPE also affect the b^c inclusive 
semileptonic decay rates, once the applied cut on the charged lepton energy becomes too 
high thus decreasing the hardness Q. These effects are conceptually similar to the biases 
in the photon energy moments, however the size of the two cannot be related to each 
other unambiguously. 

The remainder of the paper will be organized as follows: Sect. 2 presents the standard 
OPE for the photon energy spectra and their moments, based on the Wilsonian imple- 
mentation with an explicit scale separation. We address in detail the nontrivial impact of 
experimental cuts on the OPE in Sect. 3 and how they induce biases in the values of such 
truncated moments; the interplay with the perturbative corrections is studied. A step-by- 
step illustration of the computational procedure is provided. We discuss the effect of the 
non- valence four-quark operators on the inclusive B ^X s + 7 characteristics in Sect. 4. 
In Sect. 5 we comment on the corresponding bias effects in semileptonic decays before 
summarizing in Sect. 6. While the numbers and plots in the paper use the exact formulae, 
only simplified expressions are presented in the main text. Many of the full expressions 
are too cumbersome and are relegated to the Appendices. In A.l we describe in some 
technical detail the calculation of the photon spectrum in a formulation where the Wilso- 
nian prescription for the OPE has been implemented from the start. Additional non-Oy 
perturbative corrections are given in A. 2, and in A. 3 we give the numerical evaluation of 
the bias effects. 

While the present paper was finalized the preprints [14] appeared which analyzed the 
B^X s + '-f decay heavily relying on a SCET- motivated approach. We largely disagree 
with the findings of the author, which basically claim large perturbative uncertainty in 
the decays and challenge accurate calculability of the integrated characteristics derived 
from the B—*X s +i spectrum. We have added Sect. 2.1 to briefly address this controversy, 
however the detailed discussion of the SCET-based approach goes beyond the scope of 
the present paper. 
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2 Standard OPE for the photon energy moments 



The purely perturbative spectrum for b— *s+7 has been calculated in Ref. [13] (in the pole 
scheme) to first order in a s and including the BLM corrections oc (3 a 2 (f3$ = l 1— \n / = 9) . 
The results for the moments read 



(£ 7 ) pert i 



Ej>E cu t 



1 - [ai(x; fx) + <2i(x) 



hix-,11) + 6i(a;) 



(E 2 - (E y ) 2 ) pert \ El>Ecu = 
where 



a. 



[a 2 (x; fx) + a 2 {x)\ h b 2 (x; fx) + b 2 (x) 

7T L 



A) /«£ 

2 Vtt 

A) /a* x 1 
2 Vtt 



2£, 



cut 



(5) 



(6) 



Existing calculations of the perturbative corrections do not allow to specify which b quark 
mass enters x in Eq. (6), since the results in different normalization schemes differ only in 
order a 2 (without (3 ). However, on general grounds we know that the pole mass cannot 
appear in short- distance calculations; additional arguments can be found in Ref. [21]. 
Therefore, we assume mj in Eq. (6) to be the running b quark mass normalized at the 
scale fx ~ 1 GeV. 

Terms proportional to Si(x), b\(x), a 2 (x), b 2 (x) correspond to contributions from op- 
erators other than the dominant 7 . While those are noticeable for the total decay width, 
they are expected to be insignificant for the photon energy spectrum. For instance, they 
are found to shift the first and second moment by about —0.004 GeV and 0.0008 GeV 2 , 
respectively, which translates into Srrib — — 8MeV and 5 fx 2 ^0.01 GeV 2 , i.e. far below the 
accuracy of the perturbative piece alone. Thus there is no need for a detailed discussion 
of these effects. 

The values of the coefficients ai(x; fx) - b 2 (x; fx) in Eqs. (5) depend on the normalization 
point fx. They can be fixed by the fact that the total moments (including nonperturbative 
contributions) are independent of fx and by the condition that, to a given order in a s , the 
expressions at fi = reproduce the result in the pole scheme which assumes no infrared 
cutoff in the perturbative contributions. This yields 



ai(x; fx) = ai(x; 0) — -C F — C F — ^ 



4„ ^_13^ J M_ 3 
3 1 m h 18 ml 




1^ f 2 



a 2 (x; fx) = a 2 (x; 0) - ^C F ^ + ^C F —^ 



mi 



ml 



(7) 



(M is the normalization scale for a s in the MS scheme; we take it to be m b .) The analytic 
expressions for ai j2 (x; 0) and b lt2 (x; 0) are given in Appendix A.l, Eqs. (A. 15). 
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The expressions in Eqs. (7) are approximate in their /i-dependence - even to a given 
order in a s - in two respects. They account for terms only through order jm\ . Since 
n/rrib is a truly small parameter at GeV, these are tiny corrections. Secondly and 
more importantly Eq. (7) is valid only at sufficiently low E cnt . Otherwise, the simple 
prescription described above does not reproduce the actual /i-dependence of the moments: 
it introduces systematic shifts that are the perturbative analogue of the biases. In the one- 
loop approximation (including all its BLM corrections) this happens for 2E cnt >mb — 2fi, 
i.e. where the hardness Q descends below 2/i. At high E cut these 'perturbative biases' 
become significant. 

Therefore we do not rely on the simplified expressions (7), but rather compute the 
moments explicitly integrating the perturbative spectrum calculated with the Wilsonian 
cutoff at the separation scale fi. We have derived it both to order a s and to all orders in 
BLM, as described in Appendix A.l. The analytic expressions for the first-order moments 
di(x; fj,) are reasonably compact. For the BLM corrections the analytic expressions have 
also been obtained, yet they are too lengthy and we do not give them explicitly in this 
paper. We have checked that at low cuts the /^-dependence indeed follows Eqs. (7). 

The Wilsonian spectrum and, in particular its variation upon changing [i has some re- 
markable features. In contrast to the usual spectrum, it contributes even above mf,(/x)/2. 
The center of gravity of the spectrum upon varying ji does not change to the first ap- 
proximation; modifications arise in order 5fi 2 /m b , as mandated by the OPE applied to 
perturbative effects. 

The nonperturbative pieces in the energy moments are given by 

/ rp \ power 



Am 



b 




(E>-(E^r™ = § - 2p \ 2m f S + O ( — 2 ] . (9) 

We have calculated them for the local weak decay operator O7; our results agree with 
the expressions in the revised version of Ref. [7]. In principle, there are perturbative 
corrections to the coefficients of these power terms; they have not been calculated for 
realistic decays. 

The last term in Eq. (8) stands for certain effects of the nonperturbative conversion 
of a cc pair into a photon where the underlying weak decay is b-^ccs. It scales only like 
Such terms cannot arise in totally inclusive rates; they rather represent final state 
interaction among the channels mediated by O2 and O7. 

One should note that there exists no (local) OPE for such rescattering widths. Thus, 
for instance, there is no reason for m c entering there to be a short- distance charm mass. 
One may think that the meson mass Mjj is actually appropriate there; this choice is 
adopted in our numerical estimates. Fortunately, again this effect is negligible, yielding 
only 8{Ej) ~ —0.003 GeV, i.e. potentially shifting m b by as little as 6MeV. 

Without accounting for the perturbative corrections to the Wilson coefficients of 
higher-dimensional (nonperturbative) operators the total moments are given by the sum 
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of the perturbative and nonperturbative contributions. In the absence of the former, the 
energy moments obtained in the l/m& expansion do not depend on E cnt , which clearly 
cannot hold for the actual spectrum. Therefore, the thus computed moments are biased: 

(£ 7 ) ~ (£ 7 > pcrt + (£ 7 > powcr + \5m h , 
(^-(E 7 ) 2 ) ~ (^-(E 7 ) 2 ) pcrt + (^-( J E 7 ) 2 ) powcr - (10) 

Note that the biases Srrib, 5f/% are defined as the corresponding corrections to the moments 
calculated in the standard OPE in terms of the actual heavy quark parameters. We have 
introduced them with the coefficients motivated by the relation between the moments 
and the heavy quark parameters in the large-m^ limit. For this reason the biases are 
close - yet strictly speaking not identical - to the shifts to be applied to the values of m& 
and /z 2 naively extracted using the conventional prescription. In Sect. 3 we formulate the 
suggestions for evaluating these biases. 



2.1 Comments on the literature 

When the present study was completed papers [14] appeared which addressed inclusive 
b — > s + 7 decays from a different perspective motivated by the technique of the so- 
called 'Soft-collinear effective theory' (SCET). 3 Their basic statement appears rather 
paradoxical - the perturbative uncertainties were claimed to be quite significant, in conflict 
with our findings. Since we have performed explicit calculations, we conclude the large 
uncertainties in the treatment of Refs. [14] are mainly an artefact of their approach, 
intrinsically associated with the adopted SCET framework. A critical discussion of the 
latter goes beyond the scope of our paper. Instead, we shall try to elucidate the problem 
for the case at hand with minimal technical details. With SCET being a fast evolving 
field aiming at developing a uniform approach to the light-front processes in B decays 
(b — > u£is, b — > s + 7), certain controversial issues remain there. In particular, we may 
not agree with some statements in SCET. It is therefore difficult to provide an unbiased 
thorough review of the literature on the subject. As a partial solution, we refer the reader 
to the recent papers [15, 16, 14] for the novel applications and for references to the earlier 
papers and to the alternative approaches. 

The so-called Soft-Collinear Factorization in QCD is highly emphasized in Refs. [14] as 
a theoretical advantage; according to it the spectrum is represented in the form involving 
the 'hard' function H, jet function J and and the 'soft' 'shape function' S: 

— ocH-J®S, (11) 

where ® symbolically denotes a convolution. S is just the familiar nonperturbative heavy 
quark distribution function; in our Wilsonian scheme it is F(k+). In fact, the usual 

3 There are actually a number of competing versions of SCET, sometimes marked with the numbers: 
SCET-I, SCET-II, etc. A partial list of references can be found, for instance, in Refs. [14]. 



7 



factorization of the 'soft' modes employed in the Wilsonian OPE (i.e., between the effects 
from small and large distances) is sufficient for this purpose. Moreover, H ■ J is nothing 
but the (Wilsonian) perturbative spectrum d ^ cr . Contrary to the impression given in 
Refs. [14], an expression for the total spectrum in this form was presented already in 1993 
[17] (see also Ref. [18]), not to mention a number of later papers. 

Soft-collinear factorization, which goes back to the Low theorem [19] and its general- 
ization in the Gribov-Low theorem [20], in this context refers merely to how to calculate 
the perturbative spectrum itself. Regardless of anything else, it is well known how to 
calculate the perturbative corrections in B decays, and this has been done. 

The principal point of the analysis of Refs. [14] is dividing the dynamics driving inclu- 
sive b— >s+7 decays into three domains of essentially different mass scale, namely 
~ v /m f) /i hadr and ~ /ihadr- The perturbative effects residing in the first two regions were 
treated separately applying the renormalization group methods which implies a strong 
hierarchy among the scales. 

A critical discussion of the validity of a similar approach to b — > s + 7 has already 
been given in Ref. [21]. Rather than recapitulating the details here, we mention only 
a few obvious points. With the genuine Wilsonian cutoff at 1 GeV the gluon energy in 
the perturbative diagrams cannot be below 1 GeV. Kinematically, on the other hand, no 
gluon with E > ^ ~2.3 GeV can be emitted, and even those with an energy approaching 
2 GeV are strongly suppressed. This leaves a limited phase space even for single gluon 
emission and suppresses multiple emissions even more. The approximation where both 
the individual gluon momenta and their aggregate energy are much smaller than the jet 
energy m&/2, which is assumed in the Sudakov-type treatment, can never be realized in 
B decays even remotely. 

The Sudakov treatment aims at log-enhanced physical effects, it allows their resumma- 
tion. For that one employs a reorganization of the perturbative series in powers of 1/ In Q 
rather than a s , with Q denoting a supposedly large ratio of relevant scales. With a s -ha.Q 
then being of order unity, the expansion in its powers has to be resummed. Yet in B 
mesons the 'large' logs hardly reach unity. One cannot rely on an expansion valid at large 
InQ, when instead InQ < 1. The large-log series might not even converge in that domain; 
there is no reason to expect such a procedure to yield a meaningful result. This general 
problem becomes quite manifest in b^s + '-f decays [21]: for instance, the NLO terms 
dominate over the leading order turning the physically expected Sudakov suppression into 
an enhancement. 

Furthermore, the renormalization group treatment per se is inefficient in fixing the 
precise scale at which to evaluate the running strong coupling, especially in the case of 
a limited momentum interval. The concept of the renormalization invariance itself has 
nothing to say about whether, for instance, a s (A) appearing in Refs. [14] is actually 
a s (A), or a s (2A) or a s (A/2). On the other hand, a s runs fast in the considered low- 
energy domain. This is a major factor contributing to the large uncertainties of Refs. [14]. 
This ambiguity would only be resolved in the RG framework by calculating the matching 
coefficients to higher orders. At the same time, the improvement achieved directly in 
QCD by adding straightforward BLM corrections, without subdividing the perturbative 
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domain, does this far more reliably [21]. 

Let us recall that a detailed illustration of how a formally correct large-log renormali- 
zation-group treatment yields a numerically misleading result for an energy scale gap typ- 
ical in B physics, was given already 10 years ago [22] , in a similar context of the zero- recoil 
b^c transitions. That time a similar NLO renorm-group improvement accomplished by 
Neubert was claimed to yield a precision perturbative factor r]A = 0.986±0.006 [23], while 
a thoughtful application of the direct one-loop QCD calculation yielded t\a — 0.96 - a 
value later confirmed by the explicit higher-order calculations. We think the observed 
large uncertainties in the calculation of the photon spectrum in Refs. [14] root in a similar 
misconception. 

An additional advantage of our approach over SCET is that the (quite significant) l/m& 
corrections, including those in the perturbative coefficients, are automatically included; 
the KLN-type relations are respected at every order, which does not necessarily hold in 
the truncated RG approach. We emphasize that, using a commensurate language, the 
required anomalous dimensions, the BLM corrections to the 'SCET matching coefficients' 
are all automatically incorporated in the full QCD perturbative corrections we use; they 
actually can be read off from the analytic expressions of Appendix Al. 

The most accurate evaluation for B decays is obtained in the full-QCD fixed-order 
perturbative calculations. Implementing the Wilsonian approach significantly reduces 
the corrections and makes them well behaved; the effect of the running of the strong 
coupling in different kinematic domains is readily incorporated by calculating the BLM 
corrections in the same framework. This strategy had already been advocated in Ref. [21] 
and is performed in this paper. 

3 Bias corrections 
3.1 Basic ansatz 

As described in Ref. [5] and stated in the Introduction, there is no unambiguous way 
to determine precisely the biases in the photon moments, unless the nonperturbative 
dynamics for B mesons have been brought under full theoretical control. On the other 
hand with B mesons constituting the ground state we can expect the b quark distribution 
function F(k + ) describing its intrinsic motion, to be essentially a positive and smooth 
function with a pronounced peak near E 1 — rri} ) /2. There is no unique choice for such a 
function. Yet once the average and the width of the distribution are fixed by the values of 
the heavy quark parameters, the function is sufficiently constrained to allow reasonably 
accurate estimates of the cut-induced biases for naturally shaped distributions. 

For proper perspective let us recall a few basic facts. The tail of F(k+) corresponding 
to the very low E 1 and its behavior near maximal k + corresponding to E 1 approaching 
Mb/2 depend on the high- dimension expectation values and are determined by the subtle 
details of the bound-state dynamics. That implies that nothing reliable can be said about 
biases at E cut >m b /2, nor are they under full theoretical control at very low E cut . Yet for 
low E cut their size is sufficiently small to render them insignificant. Our prescription is 
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therefore meant to be applied for E cut in the range from about 1.7 GeV to 2.1 GeV. Those 
numbers follow from the expected values mf,~4.6GeV and /z^~0.45GeV 2 [2, 1]. 

The relevant parameters describing the shape of the b quark distribution function are 
its mean and variance which (in terms of E 7 ) are m&/2 and /^./12, respectively, in the 
heavy quark limit. For finite quark mass they have power corrections, which can be read 
off from Eqs. (8)-(9). For use in F(k+) we can then replace rrib and /i 2 by fn b and fj% 
with 

Zrrib bmf 
ti^l-r^&^L. (12) 

In the strict l/m b expansion one has r 3 = 1. However, the l/m b corrections are not 
negligible as is seen from the expression for the second moment. On the other hand, it 
must be positive at any mj. This is the reason to expect that the power expansion for it is 
sign-alternating, and that the l/mf correction partially offsets the last term in Eq. (12). 
Therefore, we use in practice r 3 = 0.75 ± 0.25 as an educated guess of these effects. 

Following Refs. [5, 21] we estimate the biases in Eqs. (10) using two ansatze for the 
heavy quark distribution function: 

F 1 (k + ) = Ni (A-k + ) a e ck+ 9(K-k + ), 

F 2 (k + ) = N 2 (A-k + f e~ d ^^ 6(K-k + ) , (13) 

with the parameters fixed by Eqs. (12). We will see that the two classes of functions, 
despite their different functional dependence, yield very similar results once Eqs. (12) are 
imposed fixing their mean and variance. Therefore we take their average as the central 
value, and include the difference into the estimate of the theoretical uncertainty. 

In general, the adopted distribution functions depend on the values of mj and /i 2 , and 
the biases additionally on E cut . Within the l/m& expansion the latter depend primarily 
on the combination Q = frTb — 2E cut (directly related to the hardness Q ~ Mb — 2E cut ) 
rather than on E cut and m b separately. By dimensional arguments they are actually 
functions of two variables since the absolute scale m b does not enter. Yet even with these 
simplifications the resulting expressions are cumbersome, and it would not be enlightening 
to present them here. Fortunately it is not necessary; for we have verified that the biases 
follow simple scaling behavior with reasonable accuracy sufficient for our purposes. Once 

calibrated by \J~[J%, the biases depend primarily on the ratio q = Q/yffJ^: 

S "' b =f(q,r) ~ f s (q) , = g (q, r) ~ g s (q) , (14) 

with the residual dependence on the remaining ratio r = fil/A 2 rather weak (A = Mb— rrib). 
The functions f s and g s are shown in Figs. I. (For illustration in Sect. 5 we also give 
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Figure 1: The dimensionless ratios 5mb/(/^.) 5 and \i\ as functions of the cut scale variable 
q = Q/(fi"^)2 as described in the text. The middle blue curve gives the central value, the average 
of the biases obtained using the two distribution functions. 

the corresponding estimate for the third moment of the distribution.) The difference 
between the two curves which represent the results for the two classes of functions denned 
in Eq. (13) gives a sense of the scale of the theoretical uncertainty, albeit 'cum grano 
salis': having them coincide at some point does of course not mean there is then no 
uncertainty. The functions f s (q) and g s (q) can adequately be described by the following 
simple functions for realistic values of q and r: 

f a (q) -0.415 exp [-1.3 q- 0.28 q 2 ' 5 } , g a (q) -0.80 exp [-0.60 q - 0.27 q 2 ' 5 ] . (15) 



3.2 Biases and perturbative corrections 

The above bias evaluation was based purely on the nonperturbative primordial distribu- 
tion. In reality cuts are applied to the observable spectrum resulting from the convolution 
of the perturbative and the nonperturbative pieces. Perturbative corrections modify the 
bias. Moreover, the effect of the perturbative corrections depends on their implemen- 
tation, in particular on the separation scale fi. Indeed, the moments of the combined 
spectrum are /i-independent at a given cut, as it has to be for an observable. The biases 
computed with the nonperturbative distribution alone, on the other hand, would depend 
on the choice of fi since the distribution function and, in particular m& depend on the 
normalization scale. 

It has been observed [5] that with the prototype Wilsonian perturbative spectrum of 
Ref. [21] perturbative contributions had little impact on the biases. Having calculated 
the perturbative spectrum explicitly with the Wilsonian cutoff, in this paper we adopt 
a more accurate direct approach. Namely, we determine the complete spectrum starting 
from this Wilsonian perturbative kernel. Simultaneously, we confront the results with the 
corresponding OPE predictions which include complete expressions for the perturbative 
moments derived from this spectrum, see Appendix A.l. 
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Figure 2: The perturbative effects on the biases Srrib and 5fj%. The red curve is the complete bias 
using the full spectrum. The green curve is the bias calculated without perturbative corrections, 
and the blue line shows the difference between the two curves. 

Again we find that for \i ~ 1 GeV the biases are only weakly modified by the pertur- 
bative corrections. This is illustrated by Figs. 2, where we have assumed m^{\ GeV) = 
4.6 GeV and \i\ (1 GeV) = 0.45 GeV 2 . Even so in our analysis we retain the perturbative 
shifts for the sake of completeness. At the same time, due to their smallish size we just 
approximate these shifts in brrib and in by functions of only the single variable q, 
obtained at the central values expected for the parameters: 



If necessary, these can be calculated anew at the final fitted values; the procedure con- 
verges well since the shifts are small. 

3.3 Procedure 

Let us summarize now the proposed way to obtain the numerical theory predictions for the 
photon energy moments. The inputs are the (true) values of the heavy quark parameters; 
they are assumed to be normalized at fi = 1 GeV. The procedure consists of three steps. 
First, the 'biased', i.e. uncorrected moments are evaluated by adding the nonperturbative 
and perturbative contributions. The former are given by Eqs. (8), (9). The perturbative 
ones are described in detail in Appendix A.l; they refine the previously applied Eqs. (5). 4 
The first-order integrals are given in Eqs. (A. 5) and (A. 6) with Bi denned in Eq. (A. 2). 
The similar BLM pieces of the moments, b^, are obtained by integrating Eq. (A. 11) (and, 
for b\, additionally subtracting the value of B 2 given by Eq. (A. 12)). Alternatively, one 
can use the approximate numerical interpolations Eqs. (A. 16) and (A. 17). 

Next with frit, and [i\ defined in Eq.(12) we can determine the biases Srrib and <5/z 2 
using Eq. (14): 



A pcvt f s (q) ~ 0.0103 g 3 -0.0729g 2 + 0.165g-0. Ill 
A pcrt g s (q) ~ - 0.0138 g 3 + 0.0.0437g 2 + 0.0316g-0.042 



(16) 





(17) 



The numerical difference appears insignificant, though. 
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with the dimensionless scaling variable q =Q/\J fJ% , where Q = fni ) — 2E cnt . 

Finally, we may apply additional small adjustments accounting for the perturbative 
impact on the biases. While obtained in a lengthy way using complicated expressions, 
they can be approximated quite adequately by simple polynomials as in Eqs. (16). To 
account for these effects one would replace, at the previous step 

fs(q) =► fs(q) + A^f s (q), 9.(q) =► 9.(q) + A pcrt ^(g) . (18) 

The thus corrected photon energy moments can then be confronted with the data. 

Here we illustrate the procedure using the central values of the heavy quark parameters 
[24] very close to those reported by BaBar [2] : 

m 6 (l GeV) = 4.61 GeV, /4(1 GeV) = 0.41 GeV 2 (19) 
p^(l GeV) = 0.17 GeV 3 . (20) 

Likewise we use the central number for the chromomagnetic value inferred from the B* — B 
mass difference, and we fix the one for the LS term to its theoretical expectation: 

/4(1 GeV) = 0.35 GeV 2 , p\ s = -0.15 GeV 3 . (21) 

For _E cut = 2.0GeV we have a; = 0.8677 and for the perturbative contributions we find 
from Eqs. (A.5), (A.6), (A.10)-(A.12) and (A.18): 

ai (0.8677; + ax (0.8677) = -0.0511 , &! (0.8677; + Si (0.8677) = -0.0909 
a 2 (0.8677; ^-) + a 2 (0.8677) = 0.00447 , 6 2 (0.8677; ^) + 6 2 (0.8677) = 0.00538 ; (22) 

using a s = 0.22 and (3 = 9 yields 

(S 7 ) pert U 7 > 2 .oGeV = 2.328 GeV, (E 2 - (£ 7 > 2 > per V>2.oGeV = 0.00390 GeV 2 . (23) 

The nonperturbative power corrections given by Eqs. (8) and (9) are 

(_E 7 ) powcr = -0.00715 GeV, (E 2 - (£ 7 ) 2 ) powcr = 0.0253 GeV 2 . (24) 

Adding the perturbative corrections and these power corrections results in 

= 2.321 GeV, ((E,-(E,)) 2 )\^£ V = 0.0292GeV 2 . (25) 

The last step is to correct for the biases. For our set of parameters, Eqs. (12) yield 
m b = 4.602 GeV or A = 0.677 GeV and, using r 3 = 0.75, we obtain /4 = 0.330 GeV 2 . At 

E cut = 2.0 GeV the value of the effective normalized hardness is then q = Q / \fjj^. ~ 1.05 and 

we have r = fi 2 r /A 2 ~0.72. We take the scaling form Eqs. (17) for the biases parametrized 
in Eqs. (15), and apply the shifts Eqs. (16) due to the perturbative effects according to 
Eqs. (18): 

5m b = ( f s (q) + A pcrt / S (g)) = 0.041 GeV 

5/4 = g (g s (q) + A^gM) = 0-111 GeV 2 . (26) 
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We arrive at the final bias-corrected values of the moments combining Eqs. (25) and 
(26) according to Eqs. (10): 

(£ 7 >| 2 .ogcV = (^)K + ¥ = 2.342 GeV 
<(£ 7 -<£ 7 » 2 >| 2 .0GeV = ((£ 7 -<£ 7 >) 2 >lKv - ¥ = 0.020GeV 2 . (27) 



If we follow the similar procedure for other cuts, we find, for the biased average photon 
energy and the variance, the following values for _E cut = 1.8, 1.9, 2.0 and 2.1 GeV: 

(# 7 ) IlsgcV * 2.305 GeV <(£ 7 - (£ 7 » 2 > |& ~ 0.0357 GeV 2 

(Ej) Kv - 2-313 GeV {(E 1 — {E 1 )) 2 ) |^ v ~ 0.0321 GeV 2 

(£ 7 ) IzogcV ^ 2.321 GeV ((E 7 - (E,)) 2 ) & ~ 0.0293 GeV 2 

(^ 7 ) Izi GeV — 2.329 GeV ((£ 7 - <£ 7 » 2 > & ~ 0.0271 GeV 2 . (28) 

The bias corrections read 

5m(,|i.8GeV — 0.015 GeV (fyiji.sov — 

0.0574 GeV 

5m 6 |i. 9G ev- 0.024 GeV <5/4li.9GeV - 0.0790 GeV 2 

Sm b \ 2 . 0GcV ~ 0.041 GeV 5/i 2 | 2 . 0GcV -0.111 GeV 2 

5m b | 2 .iGeV^ 0.070 GeV ^ 2 | 2 . 1GcV ~ 0.152 GeV 2 . (29) 

Then we arrive at the following final, i.e. bias-corrected predictions for the average and 
variance, 

(-E' 7 )|i.8GeV — 2.312 GeV, (( J B 7 -(E 7 )) 2 )| 1 . 8GcV ~ 0.0309 GeV 2 

1.9 Gcv — 2.325 GeV, ((£ 7 - (£ 7 » 2 > | L9Gc v ~ 0.0255 GeV 2 

(^ 7 )| 2 .ogcV- 2.342 GeV, ((£ 7 - (£ 7 » 2 > | 2 . 0Gc v - 0.020 GeV 2 

(^ 7 )| 2 .igcV^ 2.364 GeV, ((E y - (E y )) 2 ) | 2 . 1GeV ~ 0.0145 GeV 2 . (30) 

These numbers, in principle, can be compared with some published measurements (a 
more detailed analysis of the BELLE data has also been performed 5 ): 

(# 7 >|i.8GeV - 2.292 ± 0.026 ± 0.034 GeV BELLE [25] 

((£ 7 -(£ 7 )) 2 )| L8GcV ~ 0.0305 ±0.0074 ±0.0063 (GeV) 2 BELLE [25] 

(£ 7 )| 2 .ogcV - 2.346 ± 0.032 ± 0.011 GeV CLEO [26] 

((£ 7 -(£ 7 )) 2 )| 2 .oGeV^ 0.0226 ±0.0066 ±0.0020 (GeV) 2 CLEO [26] (31) 

In view of the uncertainties in the experimental numbers and their high correlations we 
refrain from detailed comparison; yet there appears a clear trend that the bias corrections 
bring the predictions into close agreement with the data, in particular for the variance. 

We want to emphasize the following features in the estimated values for biases: 
5 0. Buchmucller, private communication 
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• The biases increase steeply with the energy cuts. 

• While for E cut = 1.8 GeV they are rather insignificant, for E cut > 2.0 GeV they be- 
come comparable to or even larger than the accuracy with which the corresponding 
parameters have been be extracted from B — > X c tv. 

• The uncertainties inherent in the evaluation of biases remain rather small for m^, 
and even for ^ they are of moderate size. 



4 Effect of the four-quark operators at order OL s /m% 

As noted previously [17, 27, 21], four-quark operators (fcTjs) {sYjb) do not contribute to 
b — > s + 7 at tree level, unlike for b — > ulv. Yet hard gluon emission induces Wilson 
coefficients for them already at order a s . Here we give their contributions. Similar short- 
distance corrections affect, in general, Wilson coefficients of all other operators, including 
those determining the moments of the light-cone distribution function. The four-quark 
effects, however are usually numerically enhanced, at least where they are allowed at the 
tree level. 

To order a s the effect is expressed in terms of two expectation values 
Sv = (SikSij-lSM^- (B\b%b k s l rs j \B) 

Sa = (£|&S75& fc sty t 75S J '|£> , (32) 

where the k, I denote color indices, and yields the following contributions to the first 
three moments: 

As S T B _ x ^= 27ra s (^ + ^^4) ; (33) 

r B _ > x 7 \ m 6 3 m bJ 

the appropriate value of a s here is about 0.3. Sv,a are related to the non-valence four- 
quark parameters u>v,a, tv,a of Ref. [28] through 

S v = fccv - fry , S A = \oo A - \t a ■ (34) 

It should be noted that the contributions from the four-quark operators are l/mf- 
suppressed only in the total width. In the higher moments an enhancement arises over 
naive expectations one might infer considering the fact that the four-quark operators 
contribute to the spectrum only in the end-point region where E 1 — ^ ~/ihadr, and that 
their overall effect on the width is 1/mf. The diagrams with a gluon exchanged between 
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the light quark legs generate additional powers of t = m^ — 2E^ in the denominator in the 
forward scattering amplitude; this directly affects the higher moments, but not the total 
rate. Accordingly, we find the corrections to the rate negligible, yet not necessarily so 
for (i? 7 ) and, in particular, for the second moment determining fi^. Since contributions 
from (bTis)(sTib) can arise only via non- valence configurations in the B meson (for the 
KM favored modes), they must be somewhat suppressed; yet the degree of suppression is 
uncertain. 

The impact on the second moment and, therefore the related shift in the apparent value 
of fi^ could a priori be significant - if the generic non-valence 'color- straight' operators (in 
the terminology of Ref. [28]) are estimated through the vacuum chiral condensate {ss) 
according to 

' (B\Ws k s k (0)\B) ~ — *— (B\b l b l (0)\B) (ss) ~ -(200 MeV) 3 . (35) 



2M B 2M £ 

A dedicated discussion of the non-factorizable four-fermion expectation values can be 
found in Ref. [28]. It turns out, though that the value of Sy is expected to be particularly 
suppressed compared to the estimate (35). Sy, to leading order in 1/N C , reduces to the 
color-straight vector operator, non-valence expectation values for which (giving the s- 
quark number density at the origin) are probably small, since the density must vanish 
identically when integrated over the whole volume. The 1/N C color-octet piece of Sy 
contributing to the Darwin expectation value p^ D) is probably larger, yet still suppressed. 

A larger value is expected for the spin-dependent axial matrix element Sa, for which 
Eq. (35) can be conservatively taken as an upper bound. Yet the coefficient for Sa is 
appreciable only in (E 7 ) and in the total width where the effects are 1/™%- and l/ml- 
suppressed, respectively. Therefore we think that the related corrections in (E^) are below 
10 MeV in m b and safely at the sub-percent level for the total width. 



5 Cut-related uncertainties in b^clv 

Lower cuts imposed on the charged lepton energy in B — > X c iv decays likewise degrade 
the hardness of the transition. Yet the final state in b — > civ exhibits kinematics more 
complex than in b — > s + 7, which makes this effect less transparent. Furthermore the dy- 
namics is considerably more complicated as well, since the relevant distribution describing 
the dependence on the charged lepton energy cut, generally is neither the light-cone dis- 
tribution nor the temporal distribution describing the effects of the 'Fermi motion' in the 
small velocity (SV) kinematics. There is no reason to expect an analogous distribution to 
be universal for different moments; the power corrections are probably significant in the 
realistic settings, and probably vary depending on the moments as well. 

Nevertheless there are some common features between B^X S + 7 and B— >X c £v, in 
particular due to the conceptual similarity of the OPE treatment of both cases. Further- 
more the full impact of the degraded hardness might not reveal itself in the few lowest 
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terms in the 1/mj and a s expansion that one is, in practice, limited to in the theoretical 
description. 

As argued in Ref. [5], in this situation it is reasonable to adopt a cautious approach and 
use the bias estimates for B^X S +^ merely to assess the potential magnitude of additional 
cut-related uncertainties in the theoretical predictions for the b^civ transitions. Namely, 
the conventional OPE predictions for a particular moment, with or without the cut, 
depend on the heavy quark parameters. Assigning a cut-induced 'uncertainty' to the 
values of those parameters (even though they have definite values in actual QCD), would 
result in an additional variation in the theoretical predictions, which can be used as an 
estimate of the scale of possible cut-related exponential pieces. 

To implement this prescription in practice, we need a commensurate definition for the 
hardness in both transitions and how that scale is degraded by energy cuts. We advocate 
the following ansatz which is more accurate at low hardness: 



Q 7 ~ M B -El t -^(El t y + M|* vs. Q sl ~ M B -E^-^{E^ + M£. ; (36) 
equating the two yields the prescription 

- *&> + . . Mh - Mh .. • (37) 



for identifying equivalent energies in radiative and semileptonic decays. The resulting 
dependence is practically linear, Fig. 3. For a given Ef nt we use Eq. (37) to determine 




Figure 3: Cuts in B — > X s +~j and in B^X c tv yielding the same overall hardness Q 

the corresponding £^ ut . This allows us to estimate the size of biases associated with 
/i^, p%, ... at this hardness. Using the explicit dependence of the semileptonic moments 
in question on these heavy quark parameters, one can estimate the degree to which the 
theoretical accuracy deteriorates due to the insufficient energy release. It is interesting to 
note that Refs. [29] concluded that when the lepton cut is raised to 1.7 GeV theory looses 
control over (M|). According to Fig. 3 the lepton cut at 1.7 GeV corresponds to the cut 
on E 1 at about 2.1 GeV - and this is just the limit where the OPE for B^X s + , ~f totally 
breaks down, according to the present estimates. 
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A possible physical motivation behind this procedure has been mentioned in Ref. [5]: 
the lower cuts both on E 1 or on Eg kinematically reduce or even eliminate the con- 
tribution of high mass hadronic final states in the decay, thus introducing an effective 
'normalization scale' \l ~ Q. By virtue of the heavy quark sum rules, the integrals of 
the inclusive decay probabilities yield the values of the nonperturbative heavy quark pa- 
rameters, with the upper cutoff on the mass of the produced states playing the role of 
the normalization point for the operators [27]. For large enough values of the hardness, 
explicit perturbative corrections account for this scale dependence, yet the possible non- 
perturbative dependence at limited Q is missed. This ansatz can therefore be interpreted 
as accounting for a certain nonperturbative evolution of the heavy quark parameters at 
relatively low normalization point /i, provided it is to some extent universal, similar to 
the universality of the perturbative /i-dependence of a given nonperturbative expectation 
value at \i in the perturbative domain. 



Figure 4: The expected bias for p^ D : the dimensionless ratio ^//^/(/i^ ) 2 for the two distribution 
functions and their average (blue) 

For the sake of completeness we illustrate in Fig. 4 the expected size of the similar 
bias p 3 D in the Darwin expectation value, which following Eqs. (10) is defined (up to 
a factor —24) as the cut-related deficit in the third moment of the photon spectrum. 
As anticipated, the predictions for the bias itself are less certain, yet we see that the 
corresponding effect can be significant down to rather low cuts. 

6 Conclusions 

In the present paper we addressed in detail the theoretical evaluation of the photon 
energy moments in B — > X s +7 decays with lower cuts on the photon energy, in terms 
of the underlying heavy quark parameters. These moments are particularly important 
in connection to the precision determination of the latter. Two general aspects have 
to be addressed in calculating the moments - perturbative contributions which act as 
corrections, and the nonperturbative effects dominating the measured low moments. We 
have improved the former by computing the perturbative spectrum with the explicit 
Wilsonian infrared cutoff on the low-energy gluon modes. This has been accomplished 
analytically including the BLM corrections. 
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We believe that this approach ensures the most accurate and trustworthy results for 
the heavy flavor decays with the heavy quark mass in a few GeV range, as it is for 
the actual beauty hadrons. Extensions of the method including resummation of the Su- 
dakov radiation effects would be required for much heavier quarks, with masses exceeding 
20 GeV. 

Adding to what has been stated in Sect. 2.1, an essential element of the uncertainty 
in the numerical outcome of Refs. [14] is related to the unphysical, somewhat loose 6 
shape function S(u) arrived at in SCET in place of the truly Wilsonian F(k + ). Refs. [14] 
suggested that the properties of S(u) required for constraining its possible form, may be 
related to the true Wilson heavy quark expectation values measured, say, in B — > X c £v 
decays. Yet in practice such relations suffer from poorly controlled perturbative correc- 
tions and seem too crude. It is unclear if their accuracy can realistically be improved. 
The advantage of employing the heavy quark distribution function F(k + ) is that its mo- 
ments are directly related to the extracted expectation values; thus there is no significant 
uncertainty here. 

Accordingly, we view our calculations of the spectrum and its moments as in all re- 
spects superior theoretically to the SCET-based evaluations of Refs. [14] and as more 
accurate numerically. This conclusion is illustrated, for instance, by Fig. 5 showing the 
perturbative (Wilsonian) photon spectrum (the O7 part) in b— >s+7. It suggests that the 
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Figure 5: The significance of the higher-order BLM corrections in the Wilsonian perturbative 
spectrum (/i = 1 GeV). Blue curve is the first-order spectrum evaluated with a s = 0.35. Yellow 
curve is the usual a s plus (3oot 2 s approximation, and the (barely visible) green line is obtained 
by adding the P^a^ terms. Red curve shows the result of the full BLM resummation which can 
be performed precisely in the Wilsonian approach. 

fraction of events with Ej< 1.8 GeV is not only small, but also reliably predicted, at least 
in what concerns the perturbative effects. 

Doubts have been aired recently that the heavy-to-light transitions may not be treated 
with the conventional OPE, even abstracting from the complications associated with the 

6 The model-independent exact results for its asymptotic behavior emphasized in Refs. [14] in fact stem 
from the chosen scheme which, in some respects, is unsatisfactory itself. Thus their physical significance 
is doubtful. 
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high experimental cuts. In particular, the principal OPE results established more than a 
decade ago, were challenged. These results relate the moments of, say the photon energy 
spectrum to the local heavy quark expectation values which describe the nonperturbative 
bound-state dynamics. The OPE ensures universality - i.e. process-independence - of 
these expectation values as it should be for the characteristics of the initial bound state. 
Consequently, the same expectation values can be measured, say through the moments of 
the B^>X c lv distributions. 

We find such agnostic claims unfounded. The recent paper [4] showed that, in spite 
of the well known peculiarities of the light-like jet processes associated with the collinear 
bremsstrahlung, the perturbative corrections to the moments of the photon spectrum 
are calculable in the usual short-distance expansion, without encountering new types of 
the nonperturbative effects. Our new perturbative calculations implement this idea in 
practice. We state that the photon moments can be precisely calculated in terms of the 
analogous expectation values determined from B ^X c iv moments, as long as the lower 
experimental cut on E 1 is placed appropriately. 

On the nonperturbative side the so-called bias effects due to significant lower cuts are in 
the focus of our study. They are missed in the routinely used simplified OPE expressions, 
yet turn out to be significant in practice. We showed that for actual B decays these 
biases are affected very little by the perturbative corrections in the adequate Wilsonian 
approach with the separation scale chosen around the usual 1 GeV. This significantly 
simplifies incorporating the biases in practice. 

The main phenomenological message can be expressed rather concisely: 

• Imposing a lower cut on the photon energy in B — > X s + 7 reduces the hardness Q 
of the transition. 

• This impacts the evaluation of the energy moments in ways that are not fully de- 
scribed by the usual OPE expressions, thereby inducing systematic shifts, or biases. 

• With the presently available accuracy of the measurements those biases can no 
longer be ignored. 

• While those biases cannot be evaluated in a completely model-independent way yet, 
different ansatze for the heavy quark distribution functions, when constrained by 
the OPE, yield a consistent picture: 

- For £^ ut <1.7GeV they are more or less irrelevant. 

- For £^ ut > 2.1GeV they are in all likelihood large, yet beyond accurate theo- 
retical control. 

- In the range 1.8GeV<£2u t <2 GeV they are significant, yet can be estimated 
with reasonable accuracy. Thus they can be corrected for. 

• Lower cuts on the lepton energy in B — > X c iv conceptually have a similar impact 
on truncated lepton energy and hadronic mass moments. Yet the resulting biases in 
the values of the moments cannot be evaluated at present. We have suggested a sim- 
ple prescription to obtain semi-quantitative estimates of the potential uncertainties 
associated with the high values of E e cnt . 
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We have presented the generally rather cumbersome theoretical evaluation of the trun- 
cated photon moments in a reasonably compact form suitable even for global fits of various 
data on inclusive B decays. In our opinion, such data analyses must be done by exper- 
imentalists who can take care of numerous experimental subtleties. As the numerical 
application of our analysis, we only note that we find good agreement between the re- 
cently published BELLE data on the photon spectrum with the cut as low as 1.8 GeV, and 
the predictions based on the values of mj, /x^ and p z D extracted by BaBar from B^X c £i/ 
[2]. The data and theory predictions coincide within la of only the statistical error, cf. 
Eqs. (30) and (31), even without invoking other error bars or the uncertainties inherent 
in the theory calculations. 

Acknowledgments: We are grateful to Th. Mannel for discussions and to our exper- 
imental colleagues U. Langenegger, V. Luth and in particular to O. Buchmueller for 
invaluable exchanges. Additional comments on the difference in our treatment with the 
SCET emerged after informative discussions with M. Neubert at the Belle Workshop and 
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Appendices 

A.l Perturbative spectrum with the Wilsonian scale separation 

In Ref. [30] it has been described in detail how the required Wilsonian separation of soft 
and hard gauge field modes can be accomplished. The analysis of the semileptonic data 
based on this approach [31, 32] proved to be efficient and successful [2]. To order a s it 
corresponds to a separation based on the energy lv — \ k | of the exchanged gluon: only those 
with iv > fi are included into the Feynman integrals for the Wilson coefficients. While in 
higher orders the situation becomes considerably more involved, it can be handled for all 
BLM improvements similar to the first order case with only the relation between iv and 
\k | being modified. 

Adopting this separation scheme one computes the usual one-gluon decay probability, 
yet with the extra factor 8(\k \— jj). A massless gluon in the decay necessarily has \k\< 
hence taking /i > ^ would eliminate all bremsstrahlung. In what follows we will always 
assume 2/i < m b ; in practice, with /i around 1 GeV we have ^ ~ 0.22. The integration 
over the thus restricted phase space is straightforward and yields 

1 dF> crt (£;/i) n a s tlt _ u \ /J/-. ~ , AX „ ^ , n a 



ppcrt d £ 



Cf f £) 0(1-0 + A 5(1-0 + C F ^B 1 (£) 6'(l-£) , (A.l) 



where the second term stands for the tree-level two-body decay modified by virtual cor- 
rections. Its strength A is determined by the condition that the normalized spectrum in 
the l.h.s. integrated over the whole domain from < £ < 1, is unity. The last term with 
S'(l — £) can be viewed as the reflection of the fact that the 'kinetic' mass m b (p) in the 
presence of perturbative corrections differs from the energy of the quark at rest (a 'rest 
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energy' mass). The latter determines the perturbative end-point, in particular the photon 
energy in the two-body decay: 

4 77 rf 1 



In 



3 2 2 

The one-loop expression for the continuous part of the spectrum is 



(A.2) 



2^_ 3£ 2 -6g 1 (^-3) 



4(l-£) ^ 2(l-£) 111 v 1 w 



e < 1 - 277 



(1 - 2 ,)( 4 ^ 2 - 2; w^^ 2 ^ ) + ^ ln(2?7) e>1 _ 2?? 

These allow to readily calculate the required moments of the perturbative spectrum: 

1 1 drp crt (£;M) 



(A.3) 



Jx 



ppcrt d £ 



= < 



a;(-2a: 2 +3a;4-30) (-x 2 -2x+10) . ._, . 1 , 2 /-. x 

^ + ^ 3 Mn(l-x) + ±ln 2 (1-x) 
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for x < 1 — 2r\ 



±=|2 (-6?7a; 2 -12r/ 2 x+6?7x-6x-8 77 3 + 2 77 2 + 19 77+6) 



+ 4^^1n(2,) + 3 ln (1 _ s) + ln (2 ^ _ 1 1q 2 (2 ^ for x > 1 _ 2r/ 

(A.4) 



C F ^a 1 (x;r 1 ) = [ 

Jx 



1 1 drp crt (£;M) 

ppcrt d £ 



(1-0 ^ 



72 



(48 rf - 64 ry 3 + 72 r/ 2 - 72 r/-9 x 4 +22 x 3 +60 x 2 -96 x+23) 



+ (^(f^) ln(l-g)- Bxfa) 



^"^i^ 2 ~ (4ryx 2 +6r/ 2 x — 5r/x + 3x — 8 77 — 3) 
+ (^(f^- 8 ) In (2 V )-B 1 (r,) 



for x < 1 — 2r\ 



for x > 1 — 2r] 
(A.5) 



C F ^a 2 (x; V ) = J — 



1 drp crt (e;^) 



d£ 



C F f [tso (768 r/ 5 - 720 ?? 4 + 640 rj 3 - 480 rf + 144 x 5 - 495 x 4 - 340 x 3 

+1650 x 2 - 1020 x+61) + (1 ^ )2(3 2 a 4 2+2a ~ 17) ln (1 -x)j for x < 1 - 2?? 



^ 7T 



(l-^r) 2 (3a: 2 +2a-17) 
24 

0| (6r/x 2 +8?7 2 x — 8??x+4x — 2r/ 2 — 7r/ — 4) 



(l-~ir) 2 (l-2 7?) 



+ (l^) 2 (3, 2 + 2,-17) ln(2r?) 



for x > 1 — 2i] 
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(A.6) 



The usually considered BLM corrections can easily be calculated to arbitrary order 
with merely technical modifications. The corresponding technique has been reviewed in 
Refs. [33, 30]. First the one-gluon decay rate must be calculated with a non- vanishing 
fictitious gluon mass m g reflecting the gluon conversion into the qq or gg pairs. This 
modifies the cutoff condition: 9(\k\—fx) is replaced by 9(uo — ^) = 9(k 2 +m 2 g — /i 2 ). The 
integrals are still easily taken. All the functions additionally depend now on the new 

2 

dimensionless parameter A = (note that in Ref. [30] this ratio was denoted by A 2 ). 
Generalizing Eq. (A.l) we then set 

1 dr^(t,m g ;ri _ CF ^ mx . v)9{1H _ x)+A5{1 _ + CF ^ Bi{r) . x)5/{1 _ 



ppcrt d £ ~' n •> '/-V * V -V ~r ^ 

(A.7) 

(recall that rj< | is assumed; otherwise the kinematic constraints are modified), with 



3 



* im A) = i (1-^) -^5 + VAarcsin (f ) + ^ 



( A _2)^23T 



-arctan 



- \ (2— A) , Ail Ai f arctan 

4 V ; V ^ >y Lv / A(4-A) (r? 2 -A + l)J [VA(4-A) 

- (2 ln[r / + V / ^A]-lnA)- A2 ~^ A+2 ln( v /^-A+l + v /^A) 



7] 2 -X+l-r] . (A.8) 



v^A 



To calculate the BLM corrections, we must integrate over the gluon mass. The upper 
end of the /^-independent lower region of the spectrum now rises depending on A: the 
perturbative spectrum is independent of ji if A> 2 77 (1— £) — (1 — £) 2 : 

h (6m A)= 4(1 _ | ( ( | 2 ^ A 2 ~^ A+1) (2e 3 A 2 -3C 2 A 2 +2eA 2 +4e 4 A-14e 3 A + 17^A-lUA 
+ 4A+2e 5 -9g 4 +9e 3 + 7g 2 -15g+6) + 2gA " 2A _+f~ 3 ln (l-g ■ * A 



2(1-0 V s i-L 

at A>2r/(l-0-(l-0 2 

(A.9) 

= 4(l ( -g)"4"r ? -A) (-2g 2 A 2 +^A 2 -A 2 -4^A+12^ 2 A+4^A-9gA+A 
+8r]f + 8r] 2 f-20 7] £ 2 +4£ 2 - 12 rj 1 £+10 r? ^-8^+4r7 2 + 2 77+4) 

+ 2gA 2(l-gf 2 ~ 3 ln (2 ^ " A) at A " ^(^H W) 2 
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The above A-dependent functions replace the usual one-loop ones of Eqs. (A.1)-(A.3). 
Using Eqs. (A. 8), (A. 9) it is not difficult to obtain the BLM corrections both to the 
spectrum and to its moments. For instance, for the continuous part of the spectrum, at 
£< 1 — 2?7 we integrate the upper of the two expressions in Eq. (A. 9) over A directly from 
A = to A=l-£. At £> l — 2r7 the integral from A = to A= (l-0(£-l + 2?i) uses the 
lower expression, while the integrand in the integral from A = (1— 0(£— 1+2*7) to A = 1 — £ 
is the 77-independent upper expression. The subtraction piece employing A = 0, of course 
is integrated over A from to 00, cf. Eq. (A.12). 
Representing the spectrum as 



1 drP crt (£;/i) 



ppert 



d£ 



Uf — 

7T 



±) + m it) + - » (1-0 +A6(i-t) + 



C F - 

7T 



Br (%)+PofB 2 (%) + ...} 5>(l-0, (A.10) 



fi and Bi are given in Eqs. (A. 3) and (A. 2), respectively, while for the second-order BLM 
parts we obtain: 



38p-93Q+6£-36 6Q-31 Q+24Q-30£+18 

Ms;»7) = ^ — ^ 1^77; — 7 in(i-£) 



96(1-0 
3(0-3) 



4801-0 



16(1-0 



ln^(l-0 + ^Li 2 (0 



(A.11) 
i<l-2r, 



6r ? O-(128r/ 2 -16r/)O-(128r/ 3 -244r ? 2 +52r ? -20)£ + 52r ? 3 + 12r ? 2 -15r ? - 20 

967/(1-0 



+ 2Q-3Q-6g ln(2??+g _ 1)+ 10 ?? O-2(12r ? 2 + 3 ??+ 3)£ + 247 ? 2 -24 ??+ 6 ln(2r/) 



| 4Qr ? + 2(4r / 2 -6r ? +l)Q-2(4r / 2 -9r / + l)£-6r ? ln 2 | 
16r/0l-0 

+ fiz§ { ln (1 " ® ln (2 ^ ~ 1} " ln (1 " e) ln (2 r?) " ln (2 v) ln (2 ^ ~ 1} 

-^-w f-y )-^f-^)} «>!-*, 

The BLM correction to the coefficient of the 5'(1— can be obtained by numerically 
integrating over A 

^ = 4/ T TTM- 5l( ^ A) (a - 12) 



fi 2 



The complete /i-dependent analytic expressions for the BLM corrections to the mo- 
ments contain polylog functions at worst, but are somewhat cumbersome, and we do not 



24 



quote them here. At [i — the order-o; s spectrum and its second-order BLM correction 
reproduce the standard expressions [13] 

/^0) = ^=^ e + |E|ln(l-0 (A.13) 

/2(e;0H ; 96d-o i^w) ln(w) 

3(4 "- 3) ^ln 2 (l-0 + ?Li 2 (0) • (A.14) 



16(1-0 V v s/ 3 

Their moments, a^[x\ 0) and b^{x\ 0), respectively are given by 

, , x (-2x 2 +3x+30) (-x 2 -2x+10) , . , 1, 2/ 
/i(£0)d£=-^ ^ ^ln(l-x) + -ln 2 (l-x) 

, x -14x 3 + 41x 2 + 120x-48C(3) 1 T . , N , . 
/ 2 6 0) d£ = ^ + - Li 2 x In (1 - x 

o yfa 4 

4x 3 -31x 2 -44x-8tt 2 + 156 , ,„ , x 2 + 2x+6 T . , . 
+ gg In (1 - x) + — Li 2 (x) 

9x 2 + 18x-40+241n(x) , 2 . 4 . l T . ^ 

+ gg ^ln 2 (l-x)--ln 3 (l-x) + -Li 3 (l-x) 

(l-O^ = ^(-9x 4 + 22x 3 + 60x 2 -96x+23) + (l ~ x)(X 6 2+a: ~ 8) ln(l-x) 

/•i /A w , , ^ , 408vr 2 -369x 4 + l360x 3 -66x 2 + l380x-2305 x 3 -9x-9 T , , 
7. A(e°) 3456 + L, 2 (x) 

J 1 fi(Z;0) (l-£) 2 d£ = — ^ (l44x 5 -495x 4 -340x 3 + 1650x 2 -1020x+6l) 



(l-*)»(3*» + 2*-17) 
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f/2K;0) (i-0 2 de = 



7950 7T 2 + 7272 x 5 - 33435x 4 + 35320 x 3 - 25995 x 2 +91110 x- 74272 

86400 



-3x 4 +4x 3 + 18x 2 -36x-36 T , , (l-x) 2 (3x 2 +2x-17) , , . 
+ ^ Li 2 (x)-^ L ln 2 (1 _ x) 

(l-x)(72x 4 -528x 3 + 562x 2 -473x+1957) , . > /A , 

+ 2880 ln(1 " X) (A ' 15) 

Finally, we provide simple numerical expressions for the first and second moments of 
the second order BLM correction with [xjm^ fixed at the value 1/4.6: 



0.032-0.096y-0.67y 2 + 1.6y 3 -3.2y 4 for y > 



-0.032-0.11 y-0.33y 2 -0.47 y 3 -0.26y 4 for y < 

(A.16) 
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y = x- (1-2/4.6) 



)=/ / 2 «;A) a-o 2 de 

J a; 



0.0225-0.051y-0.066y 2 +0.040y 3 +0.24y 4 for y > 

0.0225-0.050y-0.12y 2 -0.12y 3 -0.04y 4 for y < 

(A.17) 



They are sufficiently accurate in the whole relevant domain of x between 0.65 and 0.93. 
The evaluation for different values of m b can easily be recovered from the fact that the 
above integrals are functions of the ratio /i/rrib, and at fixed the /z-dependence is well 
described by Eqs. (7). 

A. 2 N011-O7 perturbative corrections 

The non-07 spectrum was given in Ref. [13]. We have interpolated the moments of these 
contributions for ^ = 44?: 

m b 4.61 



100ai(x) = ^[0.053-0.81 (x-0.87) + 2.93 (x-0.87) 2 ] + §[-0.026 + 0.25 (x-0.87) 



+2.18 (x-0.87) 2 ] + §[0.72-8.8 (x-0.87) + 18 (x-0.87) 2 ] 
100 a 2 (x) = g| [0.0023-0.055 (x-0.87) + 0.38 (x-0.87) 2 ] + §[-0.0010 + 0.015 (x-0.87) 

-0.0413 (x-0.87) 2 ] + §[0.047-1.1 (x-0.87) + 6.5 (x-0.87) 2 ] 
100 b 2 (x) = ^[0.005-0.12 (x-0.87) + 0.71 (x-0.87) 2 ] + g[0.0007-0.042 (x-0.87) 

+0.42 (x-0.87) 2 ] + g[0.06-1.2 (x-0.87)] + 6.75 (x-0.87) 2 . (A.18) 

Since their effect is small and well below a host of the corrections left out, closer attention 
would seem superfluous. 

A. 3 Tables for biases 

Here we tabulate the values for biases obtained with the purely nonperturbative distri- 
bution functions, as referred to in Sects. 3.1 and 3.3. To make the tables compact, we 
present biases in the scaling form (14), evaluated near the expected value r = 0.75. The 
two primary entries are the averages of f s (q) (first moment) or of g s (q) (second moment) 
obtained with the two distribution functions F x and F 2 in Eqs. (13). Their interpolating 
expressions are given in Eqs. (15). To show and/or to correct for the residual dependence 
on r the third and the fourth entries, Df(q) and D g (q), respectively, linearly approximate 
this dependence of / and g: 



fs(q) — fs(q) + D f (q) (r-0.75), g s (q) — g s (q) + D g (q) (r-0.75) (A.19) 




100 61 (x) 



+0.29 (x-0.87) 2 ] + ^[0.56-7.5 (x-0.87) + 21 (x-0.87) 2 ] 
g|[0.064-0.86 (x-0.87) + 2.25 (x-0.87) 2 ] + §[0.0049-0.25 (x-0.87) 
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Table 1: The estimated values for the bias corrections around r = 0.75 



q 


fs 


9s 


Df 


D 9 


F\2 


G\2 





0.413 


0.800 


-0.0337 


0.0920 


-0.0151 


0.0209 


0.5 


0.208 


0.577 


0.00453 


0.107 


-0.00312 


0.0325 


0.75 


0.137 


0.452 


0.0141 


0.112 


0.00222 


0.0389 


1.0 


0.0854 


0.333 


0.0179 


0.112 


0.00571 


0.0440 


1.25 


0.0507 


0.231 


0.0175 


0.105 


0.00712 


0.0460 


1.5 


0.0287 


0.151 


0.0146 


0.0915 


0.00684 


0.0436 


1.75 


0.0156 


0.0938 


0.0110 


0.0743 


0.00561 


0.0375 


2.0 


0.00829 


0.0560 


0.00769 


0.0566 


0.00408 


0.0294 


2.5 


0.00225 


0.0187 


0.00329 


0.0288 


0.00168 


0.0143 


3.0 


0.000606 


0.00597 


0.00129 


0.0132 


0.000553 


0.0055 



(these dependences are approximate and were derived by varying r by ±0.25). Finally, 
the differences between the two ansatze are illustrated by the last two columns, for Srrib 
and 6(1% respectively: 

fl' 2 ^) = fs(q) ± F 12 (q) gW(q) = f s (q) ± G 12 (q) (A.20) 

where the upper and lower signs refer to the biases for F\ and F 2 , respectively. 



Table 2: The similar bias corrections for r around 1.0 



Q 


/. 


9s 


Df 


D 9 


F\2 


G\2 





0.405 


0.823 


-0.0311 


0.0760 


-0.0165 


0.0184 


0.5 


0.209 


0.603 


0.00290 


0.0932 


-0.00454 


0.0321 


0.75 


0.140 


0.480 


0.0116 


0.0988 


0.00120 


0.0399 


1.0 


0.0899 


0.361 


0.0153 


0.0999 


0.00522 


0.0465 


1.25 


0.0551 


0.257 


0.0153 


0.0951 


0.00717 


0.0501 


1.5 


0.0324 


0.174 


0.0131 


0.0845 


0.00733 


0.0492 


1.75 


0.0184 


0.112 


0.0101 


0.0701 


0.00633 


0.0440 


2.0 


0.0102 


0.0702 


0.00729 


0.0548 


0.00487 


0.0361 


2.5 


0.00307 


0.0259 


0.00332 


0.0296 


0.00225 


0.0195 


3.0 


0.000930 


0.00927 


0.00141 


0.0146 


0.000840 


0.0084 



We also give the similar table around a different value of r = (i^/A 2 = 1.0; in this case 
(r — 0.75) in Eqs. (A. 19) must be replaced by (r — 1). 
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